Strongly Interacting Two-Dimensional Dirac Fermions 
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We show how strongly interacting two-dimensional Dirac fermions can be realized with ultracold 
atoms in a two-dimensional optical square lattice with an experimentally realistic, inherent gauge 
field, which breaks time-reversal and inversion symmetries. We find remarkable phenomena in a 
temperature range around a tenth of the Fermi-temperature, accessible with present experimental 
techniques: at zero chemical potential, besides a conventional s-wave superconducting phase, un- 
conventional superconductivity with non-local bond pairing arises. In a temperature versus doping 
phase diagram, the unconventional superconducting phase exhibits a dome structure, reminiscent 
of the phase diagram for high-temperature superconductors and heavy fermions. 



Some of the most intriguing low temperature phenom- 
ena in solids occur due to strong correlations in low- 
dimensional electron systems. Prominent examples are 
the high-temperature superconductivity in the cuprates, 
the unconventional superconductivity in heavy-fermion 
compounds, or the fractional quantum Hall effect. An 
exciting novel aspect to the physics of electronic mat- 
ter arises when the electrons behave like Dirac fermions. 
The significance of Dirac-likc dispersion for understand- 
ing strongly correlated electronic systems is not well 
known. The recent breakthrough in fabricating sheets 
of graphene [l[ has provided us with a laboratory model 
of two-dimensional Dirac fermions in a hexagonal crystal 
lattice. The preparation of a well controlled, clean sys- 
tem of strongly interacting Dirac fermions in a solid state 
context appears intricate, if not impossible. In contrast 
to solid state systems, ultracold atoms in optical lattices 
can be controlled to perfection in experiments, which has 
allowed for experimental demonstrations of prototypical 
many-body phenomena, as the superfluid-Mott insulator 
transition in the Bose-Hubbard model 0, d, 0] . 

In this Letter, we show how a square optical lattice can 
be used to realize strongly interacting two-dimensional 
Dirac fermions. We employ an experimentally realis- 
tic optical square lattice [B| with an inherent effective 
staggered gauge field, recently shown to provide a Dirac- 
like single-particle spectrum at the Brillouin zone edge 
@ (see Fig. Q}t). Notably, this lattice geometry breaks 
time-reversal and inversion symmetries. We consider a 
mixture of condensed bosons and a balanced mixture of 
spin-up and spin-down fermions subjected to this lat- 
tice. The induced on-site and nearest-neighbor attractive 
interactions for the fermions are shown to be indepen- 
dently tunable, thus allowing for the study of the compe- 
tition between different strongly correlated BCS ground 
states. For temperatures around a tenth of the Fermi- 
temperature Tp, accessible with state of the art experi- 
mental cooling techniques, we find robust exotic phenom- 
ena: in addition to the conventional local-pairing s-wave 
superconductivity, an unconventional superconductivity 



Ek« k. 




FIG. 1: (color online) a) A staggered gauge field (arrows), 
induced by a time-dependent optical potential, splits the lat- 
tice into two sublattices A and B @, |(| ■ In the tight binding 
regime, the particles hop between the potential minima and 
pick up a flux ±</> around a unit plaquette. b) The single- 
particle energy spectrum in the first Brillouin zone of the 
Hamiltonian JT}. The two inequivalent Dirac points are indi- 
cated K and K'. 



with non-local bond pairing arises, which spontaneously 
breaks both, the U(l) gauge and the Cav crystal lattice 
symmetries. On the other hand, by varying the fermion 
filling fraction at a fixed coupling strength, the phase di- 
agram exhibits a superconducting dome surrounded by a 
Dirac-liquid and a Fermi-liquid at low and high doping, 
respectively. Exploring its details in a well controlled 
experimental environment may yield new insights into 
the phenomenology of high-T c superconductors or heavy 
fermion systems, which exhibit strikingly similar phase 
diagrams. 

Tailoring the Strongly Interacting Dirac Fermions: In 
contrast to other recently proposed schemes to achieve 
Dirac-like single particle spectra in optical lattices @, H| , 
we employ a square lattice with an effective staggered 
magnetic field (see Fig. [lj,) 0,0]. In the tight-binding 
limit, the kinetic energy describing tunneling in this lat- 
tice takes the form 



H 



<ij> 



h.c. 



(1) 



2 



Here, t is the hopping amplitude and 



±d> is the 



anisotropic phase factor (with the sign depending on the 
indices i and j), which denotes the presence of a dimen- 
sionless magnetic flux </> (in units of the fundamental flux 
quantum) alternating in sign across the plaquettes. The 
operators a\ (ai) and 6 J (bj) create (destroy) a particle 
in the sublattices A and B, respectively (see Fig. [I]}). 
They satisfy canonical commutation (anti-commutation) 
relations for bosonic (fermionic) particles. The spacing 
between adjacent A sites is d = A/V2, with A being the 
wavelength of the optical lattice. 

For ultracold single-component fermionic atoms, inter- 
actions are absent due to negligible s-wave scattering and 
the system is described by the kinetic Hamiltonian in Eq. 
([T]). At unit filling fraction, the low-energy band struc- 
ture possesses two inequivalent anisotropic Dirac cones at 
K = (n/d, 0) and K' = (0, n/d) along the Brillouin zone 
edges (see Fig. [T]d), with isotropy recovered for <j> = n. 
Around K, the Hamiltonian is given by 



Hq,f = 



cos 
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where a = y/2hv F , vf = s/2td/h is the Fermi velocity, 
Ok(^k) denote the Fourier-transforms of aj(bj), and the 
x and y axes are chosen parallel to the nearest neighbor 
bonds of the A sublattice. An expansion around K' yields 
a similar Hamiltonian with k x and k y interchanged. 

To implement interactions, we consider now a bal- 
anced mixture of fermionic atoms in two different sub- 
states (e.g., hyperfme states or Zeeman levels), which 
play the role of spin- 1/2 degree of freedom for ordinary 
electrons. The resulting fermionic on-site Hubbard in- 
teraction with strength Uff is tunable via the technique 
of Feshbach resonances [9(, which has been recently ex- 
ploited to realize a Mott state [n|. However, an extended 
Hubbard term cannot be easily implemented within the 
tight-binding approximation. We therefore consider an 
additional bosonic component coupling to each fermionic 
component with a strength characterized by Ubf- The 
microscopic Hamiltonian of this Bose-Fermi mixture in 
the tight-binding regime is then given by 



H — Hq f + H , 
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The first two terms arc, respectively, given by the ki- 
netic Hamiltonian (fT]) with the annihilation operators 
ai replaced by ai jCr for describing fermions with spin 
<7, or by a^B for the bosons (similarly for the sublat- 
tice B). The number operator at site i is defined by 
iij = a\ai{b\bi) for i G A(B). The Hubbard terms are 



given by H FF = (U FF /2) £\ , 



r , which denotes 



the interaction between fermionic atoms with opposite 
spins, H B b = (UBB/2)J2i n i,B(n-i,B - 1) for the in- 
teraction between bosons with Ubb > 0, and Hbf = 
Ubf^Zi „ni^ni^B for the interspecies interaction (the 



summation over sites in sublattices A and B is implicit). 
The phase diagram for mixtures of bosons with ordinary 
fermions in an optical lattice is known to comprise a rich 
phase diagram, with regimes of superfluidity, supersolid- 



ity, and phase separation 111. |12| 



To derive an effective model, we consider the ex- 
perimentally well-established regime for the bosons, 
where they form a superfluid Bose-Einstein condensate 
with condensate fraction h and healing length £ = 
(t/nUsB) 1 ^ 2 [fit- The bosonic sector is accurately de- 
scribed by the Bogoliubov theory, such that the bosonic 
density-density response function \ can be computed 



13l . Il4| . Within the parameter regime relevant for ex- 



perimentally available lattice Bose-Fermi mixtures (e.g., 
the widely used rubidium-potassium system) the consid- 
eration of the static limit of \ is sufficient [15| . To obtain 
the effective interactions for fermions, we note that the 
phonon-mediated mechanism generates, in general, non- 
local attractive interaction terms between fermions of all 
spin states, which fall off on the scale of the healing length 
£. Experiments in a bosonic 2D lattice of rubidium atoms 
by Spielman et al. Q show that typical values of £ arc on 
the order of d/y/2. With these considerations we arrive 
at the Hamiltonian H = Hq^f + Hi nt .F with 
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{U bf /U B b)xo(0 ~ Uff and g x 



(U bf /Ubb) Xi(£) are the on-site and nearest-neighbor 
interaction strengths, respectively. The static response 
functions can be readily evaluated numerically to yield 
Xo(£)i ~ 0-2, whereas the next-nearest-neighbor in- 

teractions are one order of magnitude smaller and can 
thus be neglected. Notice that the value of U bf /Ubb 
and U ff may be tuned separately, which allows for an 
independent control of both g$ and g\ . 

Competing Superconducting Phases: We now con- 
sider the superconducting (SC) instabilities that the 2D 
fermions in presence of the staggered gauge field may 
undergo as a result of the competition between the at- 
tractive on-site and the nearest-neighbor interactions 
(So 5 5i>0). While the on-site attraction favors a local 
pairing s-wave spin-singlet SC order, (a^a^) and 
(see Fig. [2b), nearest- neighbor attractions can lead to an 
inter-site bond ordering (aj|&i+ e( ,f ~ a i.lA+ei.i) involv- 
ing one of the four nearest-neighbor sites (see Fig. [5fc), 
which is called resonating valence bond [l6(. The vec- 
tors ei with I G {1,2,3,4} connect adjacent sites of the 
A and the B sublattices, i.e., e± = — e% = (x — y)d/2 and 
—e-2 = e-4 = {% + y)d/2. Within BCS mean-field theory, 
the mean-field Hamiltonian can be written as 
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FIG. 2: (color online) a) Mean-field phase diagram for the 
Hamiltonian Q with zero chemical potential, b) Real-space 
configuration of the local pairing s-wave superconducting 
phase, c) Real-space configuration of the nearest-neighbor 
spin-singlet bonding with one non-vanishing component. 
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Here, £k = 4i[cos(</>/4) cos(fc x d/2) cos(k y d/2) — 
ism((f>/4:)sm(k x d/2)sm(kyd/2)] is the full lattice 
dispersion, Eq is the condensation energy, n is 
the fermionic chemical potential, and the two 
competing SC orders are expressed as Ao,Ai ik 
in fc-space. The free energy is then given by 
F(Ao,Ai, k ) = E - fc B T£ njk m[l + exp(-£ ni k/fc B T)] 
where E n ^ with n e {1, 2, 3, 4} denote the four branches 
of the quasi-particle spectrum of the mean-field Hamil- 
tonian, which are readily obtained by diagonalizing 
2?k in Eq. ©. The variational parameters Ao and 
Ai.k are self-consistently determined by minimizing the 
free energy. We remark that the mean-field analysis 
is justified in the weak coupling regime where the 
interaction strength is smaller than the bandwidth go ft, 
gi/t < 4^2. 

The resulting phase diagram at zero chemical poten- 
tial and 7r-flux, shown in Fig. [2^,, comprises three phases: 
an unpaired phase (normal phase), an s-wave SC phase 
and a non-local pairing SC phase. The SC phases are 
separated by a first-order phase transition (dark surface 
between dashed lines) , whereas the transition to the nor- 
mal phase is second order. Due to the vanishing of the 
density of states at zero chemical potential, the Hamil- 
tonian (j4|) is quantum critical 171. Il8l|. i.e. finite critical 
coupling strengths are required to induce the respective 
BCS instabilities even at zero temperature. Note that for 
typical values of the hopping strength t ~ Er, where En 



is the recoil energy, the temperature necessary to access 
the SC phases k B T/t ~ 0(1) is on the order of 0.1 Tp, 
which is well in reach of present experiments with ultra- 
cold atoms. 

The local pairing s-wave SC phase is characterized by 
a gap function A which is constant, whereas the order 
parameter of the non-local pairing SC phase is a coherent 
mixture of orbital wave functions with a p-wave (odd) 
and a d-wave (even) symmetry, 
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The latter approximation holds around the Dirac point K 
at long wavelengths. Due to the lack of inversion symme- 
try (under interchange of the two sublattices), the order 
parameter describing the non-local pairing does not have 
a definite parity. Furthermore, the non-local pairing SC 
phase in real-space is described by a single bond formed 
with one of the four nearest-neighbor sites. This results 
in an unconventional superconductivity: not only is the 
U(l) gauge group broken in the ordered phase, but also 
the symmetry group C± v of the underlying crystal lattice. 

In the following, we consider the case of non-zero chem- 
ical potential, focusing on the unconventional SC chan- 
nel, which yields the highest transition temperature T c . 
Close to the transition, the fermion filling fraction (n) (or 
the particle/hole doping 5 = \(n) — 1|) is trivially deter- 
mined by the non-interacting limit (ignoring the Hartree 
energy) as 
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where N is the total number of sites. In Fig. [5^ we plot- 
ted T c versus g\ and S for the unconventional SC chan- 
nel. While the dispersion remains linear (Dirac-like) as 
the filling fraction (n) is slightly tuned away from unity 
(fi = 0, S = 0), the system ceases to be quantum crit- 
ical. The T c versus g\ graph develops an exponential 
tail extending all the way to g\ = (not visible for the 
scale used in Fig. [3^i) and shifts towards lower values 
of the coupling strength g\. This shift, however, does 
not grow monotonously as the filling fraction is further 
decreased but reduces again, below a filling fraction of 
approximately 0.3 (S ~ 0.7 in Fig. [3^,). Plotting tem- 
perature versus doping for a fixed value of the coupling 
strength g±, as shown in the shaded plane in Fig. [3K, 
we recognize a domc-shaped unconventional SC phase at 
intermediate filling fractions, surrounded by the normal 
phase for fillings close to unity or zero, which we termed 
Dirac-liquid (Fermi-liquid) on the left (right) side, where 
linear (quadratic) single-particle dispersion prevails. The 
dome-structure of the unconventional SC-phase is remi- 
niscent of the phase diagram for high-T c superconductors. 
The evolution of the Fermi surface, as we tune the fill- 
ing fraction from unity to zero, is illustrated in Fig. [SJs 
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FIG. 3: (color online) a) T c in the unconventional SC channel 
is plotted versus g\ for different fermion filling fractions. At 
unit filling (<5 = 0) the system is quantum critical, i.e., below a 
critical value of g\ the critical temperature is zero. Away from 
unit filling (6 > 0), T c (gi) possesses an exponentially small 
non-zero tail as gi approaches zero, which can not be seen 
in the present scale. The shaded plane at a fixed coupling 
gi/t — 1 shows a dome-like shape of the SC phase embed- 
ded in the normal phase, b) Evolution of the non-interacting 
Fermi surface as a function of doping. The shaded region is 
the filled Fermi sea. 



an unconventional superconducting phase at tempera- 
tures well in reach of present experimental techniques. 
The appearance of a dome-like structure of the uncon- 
ventional superconducting phase, when plotting the crit- 
ical temperature versus doping for a fixed value of the 
interaction, reveals an intriguing link to strongly corre- 
lated electronic materials. In fact, the staggered 7r-ffux 
Hamiltonian studied here is nothing but the mean-field 
Hamiltonian of Affleck and Marston proposed to describe 
the pseudogap phase of high-T c cuprates [22||. This in- 
dicates that our optical lattice model may play a promi- 
nent role in understanding high-T c superconductors. The 
possibility of well-controlled experimental studies of our 
optical lattice model, including the transition between 
the Dirac and the Fermi liquid, may reveal further cor- 
respondences with high-Tc cuprates or permit to identify 
missing relevant building blocks yet to be included in the 
model. 
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in the non- interacting limit. The transition between the 
Dirac-liquid and the Fermi-liquid in presence of inter- 
actions is an interesting topic for future studies, which 
might reveal further correspondences with strongly cor- 
related electron systems. The investigation of a compet- 
ing charge-density-wave instability is another important 
task deferred to further studies. 

Experimental Considerations: A promising experimen- 
tal realization of the physics described here is provided 
by a rubidium-potassium system composed of a bal- 
anced mixture of fermionic 40 K-atoms prepared in the 
\F = 9/2, m F = -7/2) and \F = 9/2, m F = -9/2) Zee- 
man components of the F = 9/2 ground state hypcrfinc 
level and bosonic 87 Rb atoms in the \F = l,mp = 1) 
ground state [l9[. The parameter U^p/UsB may be ad- 
justed via its dependence on the well depth, while an s- 
wave Feshbach resonance around 202 Gauss can be used 
to tune Uff independently. In order to experimentally 
discriminate the two SC phases discussed here, one could 
search for a signature of their distinct gap functions in 
their momentum spectra [2Cj . Correlation measurements 
similar to the one described in Ref. [2l[ should be another 
powerful method to obtain information on the nature of 
the pairing. 

Conclusions: An ultracold Bose- Fermi mixture in a 
2D square optical lattice, subjected to an effective stag- 
gered magnetic field, gives rise to strongly interacting 
Dirac fermions. A rich phase diagram arises, which ex- 
hibits a competition between a local pairing s-wave and 
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